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ABSTRACT 

The Cepheid period-luminosity relation is the primary distance indicator used 
in most determinations of the Hubble constant. The tip of the red giant branch 
(TRGB) is an alternative basis. Using the new ANU SkyMapper Telescope, we 
calibrate the TuUy Fisher relation in the I band. We find that the TRGB and 
Cepheid distance scales are consistent. 

Subject headings: galaxies: distances and redshifts - cosmology: distance scale 



1. Introduction 

The extragalactic distance scale of the Hubble Space Telescope (HST) Key Project 
( iKennicutt. Freedman &: Mould Ill995[ l is based on the Cepheid period-lumino sity (PL) re- 



lation and secondary distance indicat ors, such as t 



2000[). the supernova s tandard candle (iGibson et al 



le Tu Uy-Fisher relation (ISakai et al. 



( iFerrarese et al. Il2000[). and the funda r aental plane (IKelson et al. Il2000[). It has been crit- 



icized recently ( iTammann et al. I l2008t ISandage fc Tammann 



2000[). surface brig htness fluctuations 



2006 



20081 ) on the grounds 



that the PL relation may not be unique. Indeed, the finite width of the Cepheid instability 
strip in the HR diagram implies that nuisance parameters such as metallicity, helium abun- 
dance, and star formation history may pla y a role in determi r iing t h e PL relation . Metallicity 
was consider e d as a second parameter by iFreedman et al. I (120011 ) , ISakai et al. I (l2004j ) , and 
Macri et al. I (l2006[ l. and linear corrections were made based on measured values of [0/H] 
from the galaxies' HII regions. 



The classical extragalactic distance sca le continues to be i mportant because measure- 
ments of the Hubble Constant with WMAP (ILarson et al. Il2010l ) are model dependent: they 
integrate the scale factor all the way to high redshift. 
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It is of interest, therefore, to see how well the dist ance scale can be measured w ithout 
reference to Cepheids at all. In Papers I, II, and III (IMould &: Sakai I l2008l . l2009al Jbl) we 
calibrated the H-band TuUy-Fisher relation, surface brightness fluctuations, the Fundamental 
Plane, and Type la supernovae using the tip of the red giant branch (TRGB) ari d found a 
distance scale compatible with that of the Hq Key Project (IMould et al. 1 120001 ) . In this 
paper we use the TRGB distance indicator to examine the I-band Tully-Fisher relation. As 
discussed in Paper I, the TRGB is a good standard candle because it results from the helium 
flash on the red giant branc h, which theo r y sugg ests is relatively immune to metallicity effects 



in old stellar populations. iRizzi et al. I (120071 ) fin d hundredth of a magnit ude systematics 
in TRGB for old metal poor stellar populations. ISalaris &: Girardi I (|2005[ ) show that this 
dispersion grows substantially for young and metal rich populations. 



2. The Calibration Sample 

The calibration sample from which a Tully-Fisher relation will be constructed consists 
of 13 galaxies, which are presented in Table 1. This sample is limited to galaxies that have 
consistent methods by which their distances, I-band apparent magnitudes, and rotational 
velocities are measured. Each of these criteria is outlined below. 



2.1. Distance Moduli 



TRGB distance moduli for the calibration galaxies were taken from the Extragalactic 
Distance Database and are listed in Table 1. The procedure for obtaining distance 

moduli is described in [Jacobs et al. I ( 120091 ). CMDs for the galaxies were produced from HST 
observations, which then allows the apparent magnitude of the TRGB to be meas ured fol- 
lowing detection using the maximum likelihood algorithm of iMakarov et al. I ( 120061 ). TRGB 
distance r noduli can then be obtained using the absolute magnitude of the TRGB as cali- 
brated by lRizzi et al. I (120071 ) . and are expressed as: 



/i = mTRGB -Ai + 4.05 - 0.217[(\/ - J) - {Ay - Ai) - 1.6] 



:i^ 



for the Johnson-Cousins I and V filters, which have been corrected for galactic extinctions 
A/ and Ay. No corrections have been made for extinction within the host galaxy, but the 
observed RGB populations have been chosen such that they lie in the outer halo, where dust 



^http://edd. ifa.hawaii.edu/ 
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is scarce. For these distance moduli, la total errors are presented. 



There have been multiple measurements of the TRGB distance modulus for many of 
the individual calibration galaxies, but we adopt the EDD values due to their consistent 
calibration. 



2.2. Consistency of Distance Moduli 

Calibration of the TuUy-Fisher relation using TRGB distance moduli cannot be per- 
formed without first assessing the consistency of distance measurements to nearby galaxies 
from different methods. Firstly, we perform a two-fold consistency check between the EDD 
listed values for the TRGB distance moduli and Cepheid-based distances, as well as addi- 
tionally determined TRGB measurements. We do this by plotting in Figure 1 histograms of 
the independently determined Cepheid and TRGB distance moduli for each galaxy, along 
with the EDD measurement. It can be seen that consistency to within 0.2 mag is achieved 
for most calibrator galaxies. 



It can also be seen for NGC 3368, NGC 3621 and NGC 3627, that these /ixRGB mea- 
surements represent distinct lower estimates for their respective distances. Furthermore, the 
adopted EDD distance modulus for NGC 3351 differs from the smallest published Cepheid 
measure ment by a fu r ther .5 mag, and differs even more significantly from the /ztrgb pub- 



( 


2009 


) and 


Rizzi et al. 



( 120071 ) is the difference between I-band apparent magnitude measurement of the TRGB it- 
self, published values of which differ by ~ 1.3 mag. While this may indicate a possible 
discrepancy in the method by which the observed magnitude of the TRGB is measured, in 
both instances the ap parent magnitude of th e TRGB has been measured following maximum 
likelihood method of iMakarov et al. I ( 120061 ). We do not expect perfect agreement between 



TRGB and Cepheid distances where, as per ISakai et al. I (l2004j ). differences between these 



two measurements provides an estimate for the metallicity dependence of the Cepheid PL 
relation. Furthermore, we still expect some variation about the derived Cepheid distances, as 
each will be slightly offset from different measurements that have been calibrated differently 
(e.g. using Cepheid populations in either the LMC or the maser galaxy NGC 4258). 



In light of this uncertainty, we nevertheless choose to adopt /^trgb = 29.92 as pubhshed 
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in iRizzi et al. I (120071 ). since this exhibits a more reahstic discrepancy between Cepheid and 
TRGB-based distance measurements. 



A similar discrepancy between TRGB distance moduh mea surements exists f or the An- 
tennae galaxies (NGC 4038/39), where the values pubhshed by lSchweizer et al. I ( l2008l ) and 
Saviane et al. I (l2008[ l differ by 0.9 mag. Schweizer et al. suggest that the difference is likely 
due to the contaminating presence of AGB stars, which led Saviane et al. to misidentify 
the TRGB in the CMD. It is possible that this issue of contamination may also apply to 
measurements of the TRGB for NGC 3351. 



When considering consistency between different measurements of the TRGB distance 
modulus, /iTRGB, it is important to note that methods by which the tip itself is detected differ 
b etween groups, with re cent ly employed methods b eing the maximum likelihood approaches 



of iMendez et al. I (12002 



andlMakarov et al. I (2006), and the Gau ssian-smoothed Sob el edge 



(2001 


) and 


Sakai et al. 


(2004) 



detection filter used by 

it is beyond the scope of this paper to determine the reliability of one method over another, 
it is imperative that we adopt distance measurements that have been consistentl y calibrated. 



We ha ve favored more recent calibrations of Mj^'-'^ in light of criticisms from iRizzi et al. 



(120071 ) that a signific antly more populated CME> i s requ ired to accura tely detect the TRGB 
than those used by iLee. Freedman. &: Madore I (119931 ). And while iMadore &: Freedman 
(119951 ) suggest a minimum of ~ 100 stars for accurate detec tion, this number has been re 



vised to the far more conservative 400 - 500 stars suggested by lMadore. Mager. fc Freedman 



(120091 ). iMadore. Mager. fc Freedman I have attempted to resolve these statistical limitations 
by expressing the I-band apparent magnitude of the TRGB as a function of color (and hence 
metallicity) , such that individual stars can be corrected for c olor and essentially transferred 



to the same TRGB mag nitude. The calibration of Mj^'^^ by IMadore. Mager. fc Freedman 



agrees well with that of iRizzi et al 



Interestingly, iDalcanton et al. I (120091 ) have published multiple TRGB distance moduli 
for a number of different fields for a range of nearby galaxies as part of the ACS Nearby 
Galaxy Survey Treasury (^A NGST). Detect i ons o f the TRGB are determined using acom- 
bination of techniques from iMendez et al. I (120021 ) and lSakai. Madore. fc Freedman I ( 119961 ). 
Where the measurement of Mj^'^^ differs by less than 0.02 mag for each field within a galaxy, 
we are able to briefly assess the extent to which the choice of field within a galaxy to measure 
the TRGB magnitude affects the yielded distance modulus. We expect neither variations in 
metallicity nor internal extinction amongst halo RGB populations, where metallicities here 
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are typically low and internal extinction is assumed to be negligible. The extent of crowd- 
ing and hence contamination of the CMB by AGB stars may differ slightly between fields, 
but again we do not expect the populations of the RGB fields to differ significantly enough 
to render measurements of the TRGB inconsistent between fields. Upon inspection of the 
results for a number of galaxies, distance moduli for each galaxy vary at most by 0.1 mag, 
which are consistent to within errors for an uncertainty of 0.05 mag, which suggests that any 
variations can be largely attributed to random photometric errors. Uncertainties introduced 
by adopting the same foreground reddening v alue for each field may also contribute to the 
dispersion of distance moduli measurements ( iRizzi et al. 1 120071 ). although this suggestion 
has not been fully supported. 



2.3. Photometry 

The total sample of photometric data has been sourced from SDSS data and our own 
observations performed using the ANU SkyMapper Telescope, along with data taken from 
the literature. 



2.3.1. SkyMapper Photometry 



Photometry was obtained for four galaxies from the calibration sample (NGC 55, NGC 
247, NGC 253, a nd NGC 7793) fro m observations performed using the ANU SkyMapper 
1.35m telescope (IKeller et al. Il2007| ). located at Siding Spring Observatory, on September 
23 and 24, 2010. Observations were performed in the g and i bands and the images were 
bias-corrected and flat-fielded using bias frames and twilight fiats taken on the night. These 
were applied using the IRAF imarith package. Integration times were set at 60 s for all 
objects. 



For each galaxy, calibration from an i-band magnitud e to a standard I-band was per- 
formed using observations of standard stars from I Graham I (119821 ). 
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2.3.2. SDSS Photometry 



!^hotometry for galaxies NGC 3351 and NGC4826 was obtained from the SDSS cata- 
logue ^[_andjwas separately calibrated using the included observations of Landolt standards 
dLandolt lEoQ^ ) . 



The cahbration of SkyMapper and SDSS data is discussed in Appendix 1. Photometry 
for the remaining calibration galaxies (NGC 300, N GC 891, NGC 2403, N GC 3368, NGC 
3621, NGC 3627, and NGC 4258) was obtained from lXuUv fc Pierce I J200oh . 



2.3.3. Integrated Magnitude Profiles and Comparison of IRAF and ARCHANGEL 

Outputs 

Apparent magnitudes were measured using both the IRAF ellipse routine, and the 
ARCHANGEL photometry package developed by James Schombertlfl. ARCHANGEL per- 
forms the core routines of sky determination, ellipse fitting, and surface brigh tness profile 



generation which are common to all photometry software packages. Following ISchombert 



(120071 ). the procedure for measuring apparent magnitudes is outlined below. 



Images are first cleaned of foreground contaminating objects, and a preliminary max- 
imal radius out to which ellipses are fit is then determined. This allows for accurate sky 
background determination that avoids both foreground objects and the outer regions of the 
galaxy. In order to measure the background, boxes of a set size are placed randomly in the 
image, and a mean count is obtained for each box. The final magnitude of the background 
is then taken as the mean of each of these mean sky box values. 



Isophotal contours are then fit to the image as ellipses, which simply act as tracers 
for the stellar material. The quality of these can be manually assessed by observing rapid 
changes in parameters such as position angle and eccentricity towards the outermost ellipses. 
Once ellipses are fitted, a surface brightness profile is produced and a disk fit can be applied 
to yield apparent magnitudes. 



^ lit tp :/ /das . sdss . org/ www/html/ 

•^http: / /abyss. uoregon.edu/ ^js/ archangel/ 
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Magnitudes are evaluated at —2.5 times the logarithm of the sum of the pixel values 
within each fitted ellipse from which the total sky background inside the ellipse is subtracted. 
A final raw magnitude for the galaxy is arrived at where the integrated magnitude values 
within each successive ellipse converge, signifying that additional pixel values encompassed 
by the largest fitted ellipse are purely sky background. Thus, where these plots converge 
(see Figures 2 and 3) signifies the luminous boundary of the galaxy. 

The availability of the IRAF and ARCHANGEL fitting routines allowed us to verify the 
output of each program, particularly with regard to the integrated magnitudes. Difficulties 
arose in this particular instance of using images of nearby galaxies, where a consistent set 
of elliptical apertures could not be fit to the highly resolved (and in the case of NGC 55, 
highly asymmetrical) galactic structure. Images were either re-scaled or smoothed until a 
successful fit was achieved in either IRAF or ARCHANGEL. IRAF eUipse fits were success- 
fully apphed to all 6 calibrators for which photometry was available, and ARCHANGEL 
fits were successful on 3 of the 6. Comparisons of the fits for galaxies with both IRAF and 
ARCHANGEL fits are plotted in Figure 3. 

Agreement between the IRAF and ARCHANGEL integrated magnitudes was achieved 
for NGC 7793 to within 0.05 mag, to within 0.15 mag for NGC 4826, and to within 0.4 mag 
for NGC 55. Given that measurement errors on apparent magnitudes are typically ~ 0.05 
mag, a disagreement of 0.1 mag and above is certainly less than ideal. This disagreement in 
magnitude values can of course be attributed to the differences in the fitting algorithms that 
are employed by IRAF and ARCHANGEL. Indeed, even by visual inspection it is clear that 
ellipses are being fit out to different radii in each output. Although we have no real reason 
to doubt the validity of either fitting routine, one way to verify the output from both IRAF 
and Archangel is to visually inspect the ellipse fits and compare them to a visual estimate 
of the boundary of the galaxy. This is certainly not ideal since the entire point of eUipse 
fitting routines is to avoid such visual estimates. Nevertheless, we expect difficulty with 
these images given that most objects have their very outer edges clipped by the physical size 
of the CCD chip, along with the diffuse nature of a galaxy's outer disk, which is compounded 
by the high sky brightness under which the observations were performed. 

The difficulties encountered are particularly evident for NGC 247, whose integrated 
magnitude curve exhibits no convergence for the fitted ellipses. We can attribute this result 
to the diffuse nature of this dwarf galaxy, where the fitting algorithms struggled to define an 
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outer luminous boundary. For this particular object, we assign a larger error of 0.1 mag. 



For all remaining galaxies with a single fit in either IRAF or ARCHANGEL, we visually 
inspect the ellipse fits to ensure they have appropriately encompassed the galaxy, and assign 
a larger error estimate where there is a clear discrepancy between the fitted ellipses and a 
visual estimate of the fit. 



2.4. Photometric Corrections 



The output from photometry programs such as ARCHANGEL represents an initial, 
raw measurement for galactic apparent magnitudes. This is corrected by k times the air- 
mass where k is the atmospheric ext inction coefficient, which must be determined for each 
telescope site (ISung fc Bessell 1120001 ). Measurement errors on the observed magnitudes are 
typically given by the propagated uncertainty in the sky background measurement, which 
becomes t he dominant sourc e of error where the faint edges of the galaxy are difficult to 
determine fI Schombert 1120071 ). Sky backgrounds are determined using the skybox routine in 
ARCHANGEL, which places boxes throughout a frame, albeit away from the galaxy and 
foreground stars, and determines a mean sky value for each box. The final sky value for 
the image is then taken as the mean of each box mean, and the uncertainty taken as the 
1 / ^/n of the rms dispersion. We then apply a calibration offset to each of these measured 
magnitudes, as outlined in Appendix 1, and propagate the uncertainty in the calibration 
offset along with the aforementioned uncertainty to determine a final measurement error. 



Measured apparent magnitudes ar e then correcte d for f oreground Galactic extinction 
using the E(B— V) redenning values of ISchlegel et al. I (Il998l ). which serves as an estimate 
for the amount of light from distant objects that is scattered by dust within our own gala xy 
along its line of sight. Magnitudes were k-corrected using k/ = 0.16z, as per iHan I (119921 ). 



Magnitudes are al so corrected for internal extinction using Amj = — 7log(6/a) as per 



Giovanelli et al. I ( 1l994j ). where h/ais the axial ratio of the galaxy - the ratio of the minor and 



major axes. This internal extinction correction provides an estimate for an inclination de- 
pendent extinction, whereby light emitted within galaxies that are inclined edge-on {i ~ 90°) 
undergoes a larger amount of internal scattering relative to galaxies that are inclined face-on 
{i ~ 0°). Methods of internal extinction correctio ns are numerous, a nd do not agree as to 
how the 7 parameter should be parametrised (see ISakai et al. I l2000l for outline). Further- 
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more, some methods will in corporate a morphology - dependent offset, whereas others will 



not. The correction used by iGiovanelli et al. I ( ll994j ) is perhaps idea l then, since it ignores 



the additional morphology dependence, in light of the criticism from ISakai et al. I regarding 



the uncertainties involved in morph ology classi fication. We have chosen to evaluate the 7 



parameter as per iTuUy et al. I (119981 ) in favor of IGiovanelli et al. 



Th e final expression for the corrected apparent magnitudes is that given by IGiovanelli et al. 
Jl997ah : 

rric = rriobs - Aj + kj - Am^ (2) 



Following IGiovanelli et al. I . the total measurement error on the corrected apparent mag- 
nitudes is the quadrature sum of the measurement error and the error associated with the 
internal extinction correction. The error on the internal extinction correction is given by the 
propagation of the errors in 7 and axial ratio b/a. 



2.5. Line Widths 



For our calibration sample galaxies, we adopt the rotational velocity measurements pre- 



the pr o files taken from HI catalogues published by 



K^oriba. 



( 2005 ). lHuchtmeier et a. 
kooi ). and lKent et al. 



J2OO5I) 
j2008h . 



Theureau et a. 



Courtois et al. 



fl2006l 



■ski et al. (2004) 



Giovanelli et al. 



Springob et al. 




(2007 




Saintonge et al. 



define a window over the profile which ex- 
cludes a small portion of the integrated fiux in the wings, such that a total fiux can be 
determined within that window, and a mean fiux taken over the number of spectral chan- 
nels. The line widths are then estimated at 50% of the mean fiux density per channel. 



For galaxies with multiple profiles, a line width is measured for each available profile. 
For consistency, we adopt line widt hs derived from profiles taken from the same source as 
the cluster sample (ISpringob et al. I ), bu t this is not pos sible for our calibration sample. For 
most of the calibration sample galaxies, ICourtois et al. I line wid ths have been measured us- 
in g the either the HIPAS S BGC profiles (iKoribalski et al. II2004I ). or from profiles presented 



m 



Springob et al. I (l2005l). For NGC247 , a W^so line width was determined using the profile 



prese nted in ICarignan &: Puche I (Il990l ). For NGC 891 and NGC 2403, profiles were taken 



fromlRotsJ (Il980l 



Measurement uncertainties for each line width are limited to the spectral resolution 
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after smoothing. ICourtois et al. I provide conservative error estimates, and apply a thresh- 
old error of 20 km s~^ to indicate galaxies suitable for distance applications, in line with 



TuUv fc Pierce I fcOOOf ). 



Ot her possible line width measurements are listed in lSpringob et al. I ( l2005l ) and lGiovanelli et al. 
( 1l997a[ l. such as widths measured as 50% of the peak flux, W^q. Where consistency between 



the calibration and cluster samples is preserved, the choice of line width is largely irrelevant 
and should be motivated by breadth and accuracy of available data. We choose not to use 
the line widths included in the SFI++ sample, as we are unable to apply the same correction 
for instrumental effects to our calibration sample line widths. 



2.5.1. Line Width Corrections 



In accordance with ISakai et al. I (120001 ). line widths are corrected for inclination and 
redshift. 

= , . . (3) 

(sm«j(l + z) 

We forgo the corrections for instrumental effects and turbulence since this cannot be 
applied consistently across both data sets. Errors in the corrected line widths are given by 
the propagation of the measurement error on the observed line width and the error in the 
inclination. 



2.6. Inclinations 



One of the key parameters in Tully-Fisher applications is the axial ratio, b/a, from which 
an inclination is derived, which allows apparent magnitudes to be corrected for internal 
extinction, and to project line-of-sight rotational velocities. Following our ellipse fitting 
routines, the axial ratio is determined at the I ~ 20 mag arc sec~^ level, where our plots of 
b/a versus ellipse radius converged. Of the 6 galaxies for which w e measured inclinations , 
these values agree with values published in iPierce fc TuUv I (119921 ) , ISpringob et al. I (120071 ) 
and the RC3 catalogue (Ide Vaucouleurs et al. Ill995l ) to within 1-2°, except for NGC 55 and 
NGC 3351, where published values vary between 7° and 5° repsectively. For this study, we 
have adopted our measured values. 
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3. Fitting Procedures and the Tully-Fisher Relation 



For ou r calibration sarnple, we apply direct, inverse and bivariate least squares fits as 
outlined in iGiovanelli et al. I (Il997bl ) . Each of these approaches assumes a typical straight 
line fit for the Tully-Fisher relation of the form, M = a + blogW, where a and b are the 
zero-point and slope respectively, but each fitting method weights the residuals of the merit 
function 



X 



^ 'M,-M{\ogWf,a,by^ 



i=l 



(4) 



differently. For the direct fit, residuals in absolute magnitude are minimised such that o"j 
represents the total error in the absolute magnitude. 



For the inverse fit, line widths are expressed as a function of absolute magnitude, such 
that the dependent and independent variables have switched roles between the direct and 
inverse fits. Our model Tully-Fisher relation is now expressed as log W = M/b — a/b, and 
residuals are minimised with respect to the errors in the line width. 



For the bivariate fit, values for the slope and zero-point are determined according to 
M = a + blogW as before, but instead the merit function is now minimised according to or- 
thogonal residuals, where af = a\n + b'^a1 such that errors in absolute magnit ude and line 



width are weighted equally. Bivariate fits are performed following the method of iPress et al. 
f |l992h . 



3.1. Commentary on the Fitting Routines 



The us e of the direct, inverse, and bivariate fits are not uncommon to Tully-Fisher appli- 
cations (e.g IGiovanelli et al. Ill997bl ISakai et al. Il2000[ ). and it is certainly true that all three 
fitting procedures satisfy the basic function of providing a prediction of a variable when its 
dependent variable is known. But each regression fit will result in a different measurement 
of the slope and zero-point in cases where errors are present in both variables and are unique 
for each data point, which becomes im portant wher e Hp e xhibits particular sensitivity to 
either of these parameters. Furthermore, llsobe et al. I ( ll990l ) argue that e ach different fitting 



method represents a result with theoretically distinct implications, and iFeigelson fc Babu 
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( 119921 ) reiterate that values of Hq derived from different fitting tecliniques are not directly 
comparable. These concerns are nevertheless somewhat r nuted in applications of the TuUy- 



Fisher relation, since its origin is purely empirical (see iTuUy fc Fisher 1 119771 and sources 
therein) such that insights into the physical relation between absolute magnitude and rota- 
tional velocity remain secondary to its practicality. 



While we are purely interested in using the TuUy-Fisher relation to provide absolute 
magnitudes for galaxies with measured rotational velocities, in which instance we would be 
quite satisfied in applying the direct fit as a simple predictive tool, we must nevertheless take 
into consideration the errors in line width when fitting data, as these remain the dominant 
contribution to the overall error at smaller distances. At larger distances, the statistical 
uncertainties in the distance dominate. 



The direct and inverse fits are ap plied purely for comparison's sake to pre v ious calibra- 



tions of the TuUy-Fisher relation (e.g. lSakai et al. ll2000lJTully fc Pierce 1120001 ). Ilsobe et al. 



(Il990l ) contend that these fitting methods are applicable where the cause of the scatter about 
the fitted relation is unknown, but also note that a direct fit applicable only where the inde- 
pendent variable (rotational velocity in this instance) is measured without error. We adopt 
the bivariate fitting method as the most robust estimate of the TuUy-Fisher parameters, 
where measurement errors in both absolute magnitude and line width are treated as equal 
contributers to the overall scatter. 



3.2. A TRGB-Calibrated Tully-Fisher Relation 

Figure 4 shows the Tully-Fisher relation derived from 13 calibration galaxies, where the 
direct, inverse and bivariate fits are superimposed on the single plot, represented by the 
dot-dashed, dashed, and solid lines respectively. 

The Tully-Fisher relation for our calibration sample using the bivariate fit is thus: 

Mi = (-7.91 ± 0.62)(logVr - 2.5) - (21.11 ± 0.07) (5) 



Values for the slopes, zero-points and observed scatter for each fit are listed in Table 2. 
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As can be seen, there is a clear discrepancy between the values of the zero-point for the di- 
rect and inverse fit when compared with the value fo r the bi variate fit. But when comparing 
to the Tully-Fisher parameters of iGiovanelli et al. I ( ll997bl ) derived using the same fitting 
techniques for similarly sized samples, we can see that such disagreements are not uncommon. 



3.3. Scatter about the Tully-Fisher Relation 



We measure the scatter about the Tully-Fisher relation by determining an RMS disper- 
sion of the magnitude residuals, which corresponds to a total observed scatter, (Jobs- This 
scatter will have a contribution from the error measurements in apparent magnitude and 
distance, but what is typically f ound is that t he rn e asurement sca t ter do es not account for 
the total observed scatter (e.g., ISakai et al. I (120001 ): iMeyer et al. I ( l2008l )). In this way, an 
intrinsic scatter is implied. The intrinsic scatter is of particular interest as it illustrates an 
underlying variation between luminosity and rotational motion which can be used to explore 
how these properties are physically related. 



Following iMeyer et al. I (120081 ) . an upper limit to the intrinsic scatter can be determined 
by including in the sample to which a Tully-Fisher relation is fit only galaxies whose mea- 
surement errors are effectively minimised. From this reduced sample, the intrinsic scatter 
can be measured using, cr^^^ = <Jmeas '^^fnt- However, the lack of a robust calibration sample 
for this study means that in addition to our Tully-Fisher relation being poorly constrained, 
we have only a poor estimate of the intrinsic scatter. Since we require a measure of the 
intrinsic scatter in order to propagate this uncertainty into a final error on our Tully-Fishe r 
distance moduli, we instead adopt the intrinsic scatter of 0.25 mag via ISakai et al. I ( 120001 ). 
which was taken from their Cepheid-calibrated I-band Tully-Fisher relation. 



4. Application of the Tully-Fisher Relation to Distant Cluster Samples 

With our fitted Tully-Fisher relation, we now have the I-band absolute magnitude ex- 
pressed as a function of the logarithm of rotational velocity, such that distances to objects 
can be obtained where their rotational velocities and apparent brightnesses remain observ- 
able. 



In order to measure Hq, we require a large sample of galaxies across a range of distances 
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such that a Hubble diagram, in which an object's distance is plotted against its recessional 
velocity, can be constructed. However, we are instead interested in measuring distances 
to local and more distant clusters, rather than individual galaxies. For clusters at a large 
enough distance, we can safely neglect the physical extent of the cluster and assume that 
all galaxies assigned to a particular cluster lie at the same distance, where the foreground 
and background residuals are averaged out. Furthermore, for a cluster containing n suitable 
galaxies, we have n independent measures of the clusters distance, which allows us to reduce 
the statistical uncertainty in the cluster's distance by a factor of 1/ y/n. 

By similarly determining an average recessional velocity for each cluster, we can obtain 
a value of Hq for each cluster. The use of multiple clusters allows us to compensate for 
cluster motions relative to the Hubble flow, i.e peculiar motion, such that for a distribution 
of clusters across the sky, these motions can be averaged out. The use of multiple clusters 
additionally allows the consistency of the value of Hq to be assessed over a large range in 
distance. 



4.1. Outline of the Cluster Data Sample - SFI++ Catalogue 



The SF1++ catalogue ( ISpringob et al. 1 120071 ) . contains I-band photometry along with 
optical and 21cm line widths of ~ 5,000 galaxies. Including new photometric measurements 
for some 2000 galaxies, this catalo gue additionally consists of measurements from the SCI 
catalogue iGiovanelli et al. I (1l997a[ l. SFI, SC2, SF2, along with data published in Mathewson 
et al (1992) and Mathewson & Ford (1996). Of the total sample of 4861 galaxies for which 
reliable Tully-Fisher par ameter measurem ents were obtained, 807 galaxies across 31 clusters 
were deemed suitable by ISpringob et al. I for Tully-Fisher and Hq measurement applications. 
Included in this catalogue are measurements of the recessional velocities, Vcmb, for each 
galaxy, where these velocities have been transformed to the CMB reference frame to correct 
for motions relative to the CMB to approximate Hubble flow. 



In order to maintain consistent treatment between the calibration and cluster data sam- 
ples, we adopted frora this catalogue the galaxies whose line widths has been measured by 
Courtois et al. I (l2009l ). In emph asising this consi stency, we reject galaxies for which Wmso 



line widths were not available in ICourtois et al. I . Line widths were assigned measurement 



errors as a function of their signal to noise, where an error of 20 km s~^ or less signals a 
profile that is appropriate for Tully-Fisher applications. As a preliminary cut to the data, 
all line width profiles with an assigned error larger than 20 km s~^ were rejected. 
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We adopt all other measurements from the SFI++ catalogue where used, apart from 
the line widths for the reasons outlined above. Photometry and line width measurements for 
these galaxies are corrected in a similar manner to that of the calibration sample galaxies, 
given by equations (2) and (3) respectively. 



4.2. Selection Criteria 



Following ISakai et al. I ( |2000| ). we apply the following cuts to the cluster data: 



Galaxies with inclinations i < 40° are rejected as uncertainties in the inclination mea- 
surement begin to dominate line width errors. The deprojection of the line width to 
an edge-on measurement becomes increasingly large and uncertain as the inclination 
angle decreases. 

Both the calibration sample and the cluster should encompass the same rotational 
velocity distribution. Since the lowest rotational velocities of the calibration sample 
approach 180 km s~^, we reject galaxies from the cluster sample with smaller velocities, 
as indicated in Figure 5. 

An upper limit to the internal extinction correction of 0.75 mag is applied, correspond- 
ing to highly inclined galaxies which have more than half of their light is subject to 
self-scattering. 



Regarding the internal extinction cut-off, iSakai et al. I emphasise that the distribution 
of internal extinction values of the cluster data sample should correspond to that of the 
calibration sample. In reproducing their Figure 8 for our own data samples in Figure 6, we 
can see that an extinction cut-off to 0.75 mag remains appropriate. 



Lastly, only clusters with 5 or more galaxies and with Vcmb > 2000 km s~^ were used 
in the determination of Hq. This final sample, referred to as the Hubble Sample, consists of 
261 galaxies across 15 clusters. 



4.3. Distances and Recessional Velocities for each Cluster 



For each cluster, individual distance measurements to the constituent galaxies were 
obtained by applying the Tully-Fisher parameters to their corrected line widths to yield 
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absolute magnitudes, from which distance moduh could be determined in conjunction with 
the SFI++ apparent magnitude measurements. The cluster distance modulus is then simply 
taken as the mean of these distances, with the RMS dispersion as the corresponding uncer- 
tainty. 



From the recessional velocities of each individual galaxy, we can determine an overall 
recessional velocity for the cluster by taking the mean of these values. Whilst an uncer- 
tainty can be similarly adopted as the RMS dispersion about the mean of these values, 
we instead adopt a larger uncertainty of 300 km s~^ for each cluster's recessional velocity. 
This uncertainty corresponds to the dispersion of peculiar velocities for clusters as found by 



Giovanelli et al. I (Il998[ ). and allows us to account for any large scale peculiar motions and 



anisotropy in the distribution of clusters across the sky. 



Both of these quantities are of course dependent on the assumption that the cluster is a 
self-contained object with little substructure. This can be verified for each cluster by plot- 
ting histograms of distance moduli and recessional velocity measurements, where we expect 
these quantities to be normally distributed about some mean value. The presence of outliers 
are not likely to be indicators of measurement errors, but instead may point to errors in the 
method by which galaxies are assigned to a particular cluster. For the galaxies in th e SFI-|-- f 



template sample, cluster membership assignment follows that of iGiovanelli et al. I (Il997al ) 



in which galaxies are assigned membership based on spatial proximity to a clearly defined 
central region, and consistency with recessional velocity measurements. Galaxies that sat- 
isfy both of these criteria for a particular cluster constitute the in sample. Galaxies that do 
not satisfy spatial proximity, but have a redshift that matches that of the assigned members 



form a combined sample with the in galaxies, referred to as the in+ sample (IGiovanelli et al. 



1997aJ). For the sake of preserving as large a sample as possible, we have utilised the in+ 



sample in this study. 



Nevertheless, since we are relying on in+ sample, our cluster sample is more prone to 
outliers as we limit the sample to galaxies for which corrected magnitudes can be reliably 
determined. This becomes particularly problematic for galaxies with less than 10 galaxies, 
from which a mean distance for the cluster must be measured. However, for clusters with 
fewer galaxies, the random error measurement on the distance will be larger, and hence 
suitably reflect the loose constraint. 

With respect to measuring a mean cluster recessional velocity, we are only limited to 
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galaxies that have a rehab le cluster assignment, since the redshift for each galaxy can be 
accurately determined, independent of the properties that limit our sample by suitability of 
magnitude corrections. Again however, this method relies on the assumption that the cluster 
is a gravitationally bound system supported by the random motion of its galaxies. Under 
this assumption, the line of sight motion of each galaxy varies about some mean recessional 
motion of the cluster corresponding to random motions undergone within the cluster. Any 
evidence of substructure, in which distinct groups that have yet to coalesce would have their 
own mean cz values, would be revealed by a cz histogram with multiple distinct peaks. 



For this study, recessional velocity measurements for each cluster have been determined 
us ing the entire SFI++ temp late sample, whic h were t hen compared with values published 



m 



Masters et al. I (120061 ) and iGiovanelli et al. I (Il997bl ). This is to ensure that the sample 



from which we are working is suitably large enough for accurate measurements to be made. 
The immediate limitation of course in using only these galaxies is that we have only sam- 
pled the spiral galaxies in each cluster, and ignored the large number of elliptical galaxies. 
We may be able to take some solace in the fact that since we are sampling spiral galaxies 
only, which are found in greater numbers in the outer regions of most clusters, we avoid the 
large infall motions of galax ies closer to cluster cores which can induce large dispersions in 
an overall cz measurement (iMould et al. II2OOOI ). All our measurements are consistent with 
pubhshed values, apart from cluster A2634, where we instead adopt the published value in 
favo ur of our own. F or A2634, our measurement is based on 22 galaxies, whereas the value 
from [Masters et al. I is based on 200 galaxies, thereby providing a better sample for analysis. 



Final distances, recessional velocities, and values for Hq for each cluster are presented 
in Table 3. For each parameter, la random errors are included. 



5. Hubble Diagram and the Value of the Hubble Constant 

With distance and recessional velocity measurements, Hq values can be determined for 
each individual cluster, and a weighted mean can be adopted, where the random errors are 
used as the weights. For the data sample presented in Figure 7, the value of Hq arrived at 
via a TRGB-calibrated TuUy-Fisher relation is 79 ± 2 (random) km s~^ Mpc~^. 

We can further construct a Hubble Diagram (Figure 8). 
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5.1. Error Propagation and a Final Value for Hq 

Random and systematic errors are propagated independently into two final uncertainties 
for Hq. Random errors can be reduced statistically with relative ease, but systematic errors 
cannot be reduced in this manner, and hence require refinement of existing methods, or new 
methods altogether. These errors must be tracked carefully through each stage of the cal- 
culation from the calibration sample magnitudes and line widths until a value of Hq is derived. 



When determining a total error for a particular quantity in which we are propagating 
multiple individual errors, these are typically added in quadrature since these measurement 
uncertainties are independent of one another. The linear addition of errors instead indicates 
a correlation between errors, such that an overestimate of one quantity implies an overesti- 
mate in the other. But where error measurements are independent, quadr ature addition of 
uncertainties provides a smaller, and more appropriate overall uncertainty ( iTaylor Ill997l ). 



The following error propagation recipe follows closely that of ISakai et al. I (j2000f ). 



5.1.1. Errors in the Calibration Sale 



Random errors in the calibration scale are propagated as an uncertainty in Hq through 
the intrinsic dispersion of the Tully-Fisher relation. Systematic errors arise in the calibration 
scale through the zero-point error in MT^*^^, and th e zero-point uncertainty in the Tully- 
Fisher relation. As determined by lRizzi et al. I ( 120071 ). the systematic uncertainty in Mf^*^^ 
is 0.02 mag. 



5.1.2. Random Errors in the Tully-Fisher Parameters 

Random errors in the corrected apparent magnitude are given by the quadrature addi- 
tion of the uncertainty in the measurement and the internal extinction cor rection, the latter 



of whic h propagates the errors in the ellipticity and 7 parameters. As per iGiovanelli et al. 



(1l997al ). we assume an uncertainty of 25% for 7. Uncertainties in the redshift correction and 
galactic extinction are sufficiently small so as to be ignored. The final expression for the 
error in mc is given by: 

d = + (0.251og(l - e))' + (l^^^e) (6) 
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Errors in the corrected line widths are purely random, and are given by the quadrature 
addition of the uncertainty in the measurement and the inclination: 



-w 



1 



'1 + z)smi 



[1 + z){smiy 



(7) 



where the uncertainty in the inclination, ej is the propagated uncertainty in the ellip- 
ticity: 



1 -e 



(i-e)^-gg 



The uncertainty in the ellipticity is given by the era p irically determined relation, eg = 
0.09 — 0.12e + 0.037e^, as presented in iGiovanelli et al. I (Il997b| ). It should be noted that 
this evaluation of eg has the potential to underestimate the error for low ellipticities, but 
for the most part these are not present in the sample, since low ellipticities translate to 
approximately face-on galaxies which are excluded from the sample. Thus, a final error on 
logW is given by: 

^logW — ri ^o/iTJ^ 



5.1.3. Random Errors in the Tully- Fisher Distance Moduli and Hq 

For each individual galaxy, the random error in the distance modulus is the quadrature 
addition of the errors in the apparent magnitude measurement and the random error in the 
absolute magnitude. The random error in the absolute magnitude consists of the uncertain- 
ties propagated by the Tully-Fisher relation: the intrinsic dispersion, cxint, and the line width 
error, where the error in logW is multiplied by the slope of the Tully-Fisher relation, b, to 
be expressed in magnitudes. 

^1.,rand = + ^^logW + ^int (10) 

For each cluster, a mean distance modulus is adopted as the cluster distance, for which the 
uncertainty is simply taken as the RMS dispersion. There is a 1/ ^Jn dependence on this 
error, and hence this random error component is reducible for clusters with a larger number 
of galaxies. 
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For individual cluster measurements, the random error in Hq is comprised of the errors 
in the distance scale and the error in the adopted recessional velocity, ecz, of the cluster. The 
uncertainty on the measured cz is simply the dispersion of the individual member galaxy cz 
values, divided by the square root of the number of galaxies in the cluster. However, this 
uncertainty reflects only the extent to w hich the averaged mot ion of the galaxies represents 



motion of the cluster as a whole. As per [Masters et al. I (120061 ) . we adopt an uncertainty of 



300 km s ^ for the recessional velocity of e ach cluster to ac c ount for the peculiar motion 



of each cluster, which is in accordance with iGiovanelli et al. I (119981 ). The error in Hq for a 
particular cluster is then given by: 

= (^)' + iOA6Roe,,randf (11) 

where d is the distance to the cluster in Mpc, determined from the average distance modulus. 

We adopt as a final random error for Hq as the RMS dispersion about the weighted 
mean, reduced by a factor of l/y/n, where we have used n = 15 individual measurements of 
Ho taken from 15 clusters. 

5.1.4- Systematic Uncertainty in Hq 

The systematic error in the calibration scale, consisting of the zero-point error on Mj^'^^ 
and the zero point error in the Tully-Fisher relation propagate in to a final systematic 
uncertainty in Hq as a systematic distance error. This is defined as: 

^H,syst = 0.46Hoe^^sj/st (12) 

5.2. A Final Measurement and Uncertainty for Hq 

Following the above procedures, we determine our final value for Hq: 

Ho = 79 ± 2 (random) ± 3 (syst.) km s"^ Mpc"^ 
5.3. The Sensitivity of Ho to Changes in the Tully-Fisher Parameters 



Our least squares regression fits to our calibration sample are of course susceptible to the 
small number of calibrating points, where for such small samples outliers have a significant 
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impact on the fit. As such, were this cahbration sample to be extended, we would expect 
the slope and zero-point to change non-negligibly. Therefore, it is important to assess how 
changes in these Tully-Fisher parameters affect our final value for Hq. Using the values of 
the slope and zero-point for the bivariate fit, we evaluate the change in Hq for a Sa change in 
eacl]0, in order to anticipate possible values of Hq that could be achieved with a more robust 
fit. For a ± 3cr change in the slope, our value of Hq changes by 2%, and for a ± Scr in the 
zero-point, Hq changes by 9%. As can be seen by this quick analysis, the final value of Hq 
exhibits minimal sensitivity to the slope, but a much stronger sensitivity to the zero-point. 



6. Conclusion 

In this paper, we have utilised the TRGB standard candle to construct a distance 
scale by which the expansion rate of the universe, or the Hubble constant, is measured. 
By determining apparent magnitudes for 11 galaxies with known TRGB distance moduli 
using observations from the SDSS catalogue, the ANU SkyMapper telescope along with data 
adopted from the literature, we have constructed a Cepheid-independent I-band Tully-Fisher 
relation: 

Mi = (-7.81 ± 0.65)(logiy - 2.5) - (21.11 ± 0.07). (13) 
This Tully-Fisher relat ion was then applied to the SFI++ Tully-Fisher template data 



from lSpringob et al. I ( 120071 ). which initially consisted of 807 galaxies across 31 clusters. After 
applying cuts for line width and line width measurement error, internal extinction, inclina- 
tion, and a minimum cluster number, we maintained a sample of 261 galaxies across 15 
clusters. Applying the Tully-Fisher relation to line widths of these galaxies yielded absolute 
magnitudes for each galaxy which, when combined with their listed apparent magnitudes, 
gave a distance to each galaxy. For each cluster, an average distance could then be obtained, 
in addition to an average recessional velocity from measurements of individual galaxy red- 
shifts. These two measurements provide the basis for measuring the Hubble constant via 
Hubble's Law, cz = Hor. 



Random and systematic errors have been propagated through each step of the distance 
scale into a final uncertainty for Hq. Systematic errors continue to dominate distance-scale 
based measurements of Hq, which are principally introduced in the calibration scale, al- 
though we have reduced the systematic error effectively using the well constrained absolute 



^The 3cr do not represent confidence intervals, but simply possible variations. 
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magnitude of the TRGB in favour of the Cepheid distance scale. The final value of Hq is 
taken to be 79 ± 2 (rand.) ± 3 (syst.) km Mpc~^, which represents an uncertainty 
of 4%. While our measurement of Hq repre sents a upper esti mate relative to other recent 
determinations (74.2 ± 3.6 km s"^ Mpc"^ jRiess et al. II2009I ): 73 ± 2 (ra nd.) ± 4 (syst.) 



kms-i Mpc"^ (IFreedman fc Madore 1120101 1: 73 ± 5 (rand.) km s"^ Mpc"^ (IMould fc Sakai 
20081 )). but it is still consistent given the uncertainties. 



Agreement with the results of iFreedman et al. I (120011 ) implies that the TRGB and 
Cepheid distance scales are consistent. 
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This research has made use of IRAF, which is distributed by NOAO. NOAO is operated by 
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Appendix: Calibration of Photometric Data 



As outlined in lBessell I (|2005[ ). the idea behind calibration is to simply place photometric 
measurements from multiple instruments onto a single scale. 



A. SkyMapper Images 



In calibrating the SkyMapper images, standards from the El region from I Graham 



(Il982l ) were observed only once on September 24, 2010. Three exposures each were taken 
of the field in SkyMapper's g and i bands, from which an average instrumental magnitude 
for each star was determined. Magnitudes were determined using the IRAF phot pack- 
age, and corrected for atmospheric extinc tion using the k-coefficients for the Siding Springs 
Observatory from ISung fc Bessell I (120001 ) . These were science verification observations for 
SkyMapper hence the limited calibration available. 



Despite the small sample of standards, we have obtained a calibration applicable to a 
suitable range in color. Fitted to the data are least squares regressions, one minimising resid- 
uals in the offset, and the other minimising residuals in both the offset and the color. Errors 
in the color are the propagated errors from the individual magnitude measurements, where 
the errors in the instrumental magnitude are simply the RMS dispersion about the adopted 
mean since multiple measurements were available. Errors in the offset measurement are the 
quadrature addition of the errors in the I-band magnitudes from Graham and the errors on i. 

Also superimposed on the plot is the weighted mean which, over the plotted range in 
color, is consistent with the least squares regressions to within 0.05 mag. Added to the 
fact that there is no clear trend with color, we adopt a constant offset allowing conversion 
between SkyMapper's instrumental i magnitude and the standard I magnitude. Thus, we 
determine an overall calibration offset constant of - 3.31 ± 0.05 mag, where we have taken 
the RMS dispersion as the error. 



B. SDSS Images 



Standard fields for the SDSS cat alogue we r e obta ined through the SDSS Data Archive 
Serveil^l and I-band magnitudes from Landolt I fll992l lR were used. Atmospheric extinction 
coefficients for the Apache Point O bservatory (where the SDSS observations were performed) 
are taken from iHogg et al. I ( 120011 ) . 



In an identical manner to the SkyMapper calibration, Sloan instrumental colors g — i 
for each star were plotted against an offset, I — i, in Figure 10. Superimposed on this plot 
are the weighted mean along with a least squares regression. 

As this figure indicates, the mean is consistent with the least squares fits to 0.005 mag 
over the covered color range, such that we can suitably apply a mean offset of 2.39 ± 0.01 
mag to calibrate all SDSS images. 



^http://das. sdss.org 

^http://www.cflLt. hawaii.edu/ObsInfo/Standards/Landolt/ 
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Table 1. Calibration sample of galaxies 
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29.08 ± 0.06 


8.43 ± 0.05 %o 


3627 


11 20 15.0 


+12 59 30 


0.002425 


0.032 


3 


0.39 " 


345 ± 15 


29.59 ± 0.09 


7.71 ± 0.05 "1 


4258 


12 18 57.5 


+47 18 14 


0.001494 


0.016 


4 


0.64 " 


385 ± 20 


29.41 ± 0.04 


7.22 ± 0.05 " 


4826 


12 56 43.69 


+21 40 57.5 


0.001361 


0.04 


2 


0.48 


319 ± 16 


28.2 ± 0.03 


7.21 ± 0.03 


7793 


23 57 49.83 


-32 35 27.7 


0.000757 


0.018 


7 


0.34 


177 ± 13 


27.79 ± 0.08 


7.91 ± 0.07 



Note. — Presentation of the parameters for our calibration sample of galaxies. Values 
for ellipticity and measured I-band niagnitude are fr o m th is paper unless where specified. 
Values marked with " are taken from iTuUy &: Pierce I ( 120001 ). 



Note. — In column (1), the galaxy's NGC number is listed. In columns (2) and (3) are the 
right ascension and declination of the object in J2000 coordinates. The redshift, foreground 



reddening value (taken from ISchlegel et al. I ) , and the morphological T code are listed in 
columns (4), (5), and (6). Ellipticity values are listed in column (7), with the co rresponding 
source code. The measured rotational velocities adopted from ICourtois et al. I (apart from 
NGC 247, NGC 891 and NGC 2403) are listed in column (8). Column (9) contains the 
TRGB moduli as taken from the EDD (apart from NGC 3351). Column (10) contains the 
measured I-band magnitudes, with the corresponding source code. 
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Table 2. Tully-Fisher pcirameters. 



Fit 


b 


a 


(^obs 








(mag) 


Direct 


-7.37 ± 0.26 


-21.02 ± 0.03 


0.54 


Inverse 


-9.97 ± 0.67 


-21.04 ± 0.15 


0.63 


Bivariate 


-7.91 ± 0.62 


-21.11 ± 0.07 


0.52 
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Table 3. Cluster data 



Cluster 


N 






1^ 








D 


VcMB 


Vcmb/D 








(mag 


.) 




(Mpc) 


(km s~^) 


[km s Mpc ) 


N383 


21 


33, 


.95 


± 


0, 


.11 


62 


± 3 


4718 


77 ± 6 


N507 


13 


33, 


.90 


± 


0, 


,13 


60 


± 4 


4660 


77 ± 7 


A194 


7 


33 


.98 


± 


0, 


,13 


62 


± 4 


5002 


80 ± 7 


A262 


32 


34, 


.01 


± 


0. 


,09 


63 


± 3 


4710 


74 ± 6 


A397 


5 


35, 


.15 


± 


0. 


,40 


107 


± 20 


9746 


91 ± 17 


A400 


41 


34, 


.63 


± 


0. 


,05 


85 


± 2 


7011 


83 ± 4 


Cancer 


35 


33 


.94 


± 


0. 


,07 


61 


± 2 


4976 


81 ± 6 


A779 


12 


34, 


.98 


± 


0, 


,10 


99 


± 4 


7255 


73 d= 4 


Antlia 


5 


32 


.89 


± 


0, 


,28 


38 


± 5 


3142 


83 ± 13 


Hydra 


9 


33 


.59 


± 


0, 


,10 


52 


± 2 


4118 


79 ± 7 


A1367 


16 


34, 


.65 


± 


0, 


,10 


85 


± 4 


6711 


79 ± 5 


Coma 


27 


34, 


.40 


± 


0. 


,08 


76 


± 3 


7070 


93 ± 5 


A3574 


5 


33 


.83 


± 


0. 


,29 


58 


± 8 


4869 


83 ± 12 


Pegasus 


20 


33 


.49 


± 


0. 


,10 


50 


± 2 


3447 


69 ± 7 


A2634 


13 


35, 


.18 


± 


0. 


,09 


108 ± 5 


8895 


82 ± 4 



Note. — Presentation of the final Hubble sample cluster data. For each cluster, its average distance 
modulus and corresponding distance in Mpc, along with CMB recessional velocity and the resultant value 
of Ho are listed. 
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Fig. 1. — Distribution of distance moduli measurements for each of our calibrator galaxies. 
Cepheid-based measurements are plotted in blue, TRGB-based measurements are plotted 
in red, and our adopted EDD TRGB values are plotted in black. Width corresponds to 
measurement error, and vertical scaling is arbitrary. 
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Fig. 2. — Calibrated integrated magnitudes are plotted as a function of isopliotal (semi- 
major) radius for each galaxy. The dashed line shows the final adopted magnitude. Fits for 
each galaxy were performed in IRAF. 
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Fig. 3. — Integrated magnitudes determined using IRAF are plotted on the left, with 
ARCHANGEL magnitudes on the right. The dashed line indicates the adopted raw mag- 
nitude. For NGC 7793, we have additionally plotted the integrated magnitudes determined 
using the ARCHANGEL extremelsb fitting routine in red, which accommodates faint galax- 
ies. 
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Fig. 4. — Measured I-band TuUy-Fisher Relation. 
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Fig. 5. — Distribution of corrected line widths for the cluster sample. The vertical dashed 
line marks 180 km s~^. 



- 38 - 



0.6 0.8 
Internal Extinction (mag) 



Calibration Sample 
Cluster Sample 




1.4 



Fig. 6. — Comparison of internal extinction distributions for the calibration and cluster 
samples. 
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Fig. 7. — Values for the Hubble constant for each individual cluster (with la random errors), 
where the weighted mean is plotted as the dashed line. These values are published in Table 
3. 
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Fig. 8. — Hubble diagram for the sample of clusters. The slope of the dashed line corresponds 
to the weighted mean for Hq. 
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9. — SkyMapper I-band calibration plot. 
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Fig. 10. — SDSS I-band calibration plot. 



